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A novel D-model of wave turbulence is presented which allows to reproduce in a single frame 
various nonlinear wave phenomena such as intermittency, formation and direction of energy cascades, 
possible growth of nonlinearity due to direct energy cascades, etc. depending on the initial state. No 
statistical assumptions are used, all effects are due to the behavior of distinct modes. Classical energy 
spectra E ~ 0~ v , v = const > 0, for dispersion function of the form ~ k a , a > are obtained 
as a particular case of a more general form of energy spectra: E ~ 0~" , 2 + a -1 < v < 2(2 + a -1 ) 
where magnitude of v is defined by the parameters of the initial excitation. D-model is a generic 
model which can be expanded into a hierarchy of more refined models including dissipation, forcing, 
etc. D-model can be applied to the experimental and theoretical study of numerous wave turbulent 
systems appearing in hydrodynamics, nonlinear optics, electrodynamics, convection theory, etc. 



1. Introduction. 

Stationary energy spectra in the statistical theory of wave 
turbulence follow the celebrated Kolomogorov-Zakharov 
(KZ) law Ek ~ k~ u , with dispersion relation oj ~ k a , 
v, a > being constants, and k the wave number. These 
spectra are obtained under certain assumptions - like 
the smallness of non-linearity, infinite boxes, existence of 
an inertial range (k\ < k < k 2 ), local interaction in k- 
space etc. - in close analogy to ordinary turbulence by 
replacing vortices with waves. 

Attempts to verify this theory led, however, to contro- 
versial results (for a review see [2|), among them in sev- 
eral cases the lack of the generation of the energy cascade 
and instead the production of regular wave patterns, 
Moreover, if a cascade occurs, its spectrum consists of 
two distinct parts: a discrete spectrum and a continuous 
spectrum. The former, and frequently even the contin- 
uous spectra do not follow the KZ law. Well-known ex- 
amples are the elastic thin vibrating steel plate, Q, and 
gravity waves in mercury, [f|. Also, surface waver waves 
produced in laboratory in a flume of size 12 x 6 x 1.5 m 
developed a strongly non-linear discrete energy cascade 
lacking any continuous spectrum, [f|. The form of the 
energy spectra depends, in addition, on the parameters 
of the initial excitation, 0, [|| . In experiments with capil- 
lary waves in He strong non-locality was found; the local 
wave-amplitude maximum appeared at frequencies u> of 
the order of the viscous cut-off, . 

Attempts on improving the statistical wave turbu- 
lence theory either refer to frozen [lpj . s andp ile fTTj ) . 
mesoscopic fl2T | or finite-dimensional models [13| . neither 
allowing for the inclusion of the dependence on the 
initial conditions or finite sizes. In the present Letter 
we propose a novel model (called further on D-model) 
based on discrete wave turbulence 0, EH and including 
energy cascades. It provides a uniform frame for study- 
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FIG. 1: Color online. Nonlocal interactions in a 4- wave sys- 
tem, ~ k a . Each couple of (red) dashed lines with equal 
lengths correspond to specific choice of a parameter t. 

ing finite-size effects and deducing the generic form of 
turbulent energy spectra depending on the parameters 
of initial state. 

2. D-model. 
2.1. Excitation of a (quasi-) resonant mode. 

2.1.1. Within the discrete model of wave turbulence 
[16] the solutions of the resonance conditions form a 
group of independent resonance clusters, each cluster be- 
ing graphically represented by its non-linear resonance 
diagram (NR-diagram) which explicitly defines the dy- 
namical system of the cluster and its conservation laws. 
The solutions of the corresponding dynamical system 
may be either regular or chaotic, depending on the struc- 
ture of the cluster and the details of the initial excitation. 
This holds for quasi-resonance with small frequency mis- 
matches. 

2.1.2. (Quasi-) resonant modes occur all-over k-space; 
they are not restricted to the range k < k\. 

2.1.3. (Quasi-) resonance interactions may not be lo- 
calized in fc-space. In a 4-wave system with dispersion 
function ~ k a modes with arbitrary large wavelength 
difference can interact directly (see FigJTJ). In this case 
a parametric series of solutions of resonance conditions 
kf + k 2 = &3 + ki + k 2 = &3 + &4, can be eas- 
ily written out as k\ = (k x ,k y ), k 2 = (t,—k y ), k$ — 
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(k x , —ky), k^ = (t, k y ), with real parameter t. 
2.1.4. In addition, non-resonant modes, i.e. 



modes 
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of large frequency mismatch, do exist. In fact, most of 
the modes in a 3-wave system are indeed non- resonant. 
These modes do not change their energies on the corre- 
sponding resonant time-scale. 
2.2. Excitation of a non-resonant mode. 

In this case the generic mechanism of instability can 
be described as an interaction of three monochromatic 
wave trains 0+ + 0_ = 20o where 

0+ = 0o + A0, 0_ = O - A0, < A0 < 1. (1) 

This type of instability is quite general and is known in 
various areas of physics under different names, e.g. para- 
metric instability in classical mechanics, Suhl instabil- 
ity of spin waves, Oraevsky-Sagdeev decay instability of 
plasma waves, modulation instability in nonlinear optics, 
Benjamin- Feir instability in deep water, etc., (l 61 ] . 

Conditions for the modulation instability to occur have 
been written in terms of the increment of instability J; 
its form may differ for diverse wave systems and for the 
same wave system with different magnitudes of the non- 
linearity parameter (e.g. [l7| and [l8j . for surface wa- 
ter waves). However, as I can always be presented as a 
polynomial on 0, k and A0, it does not affect our general 
scheme (presented below) for computing energy spectra. 
To shorten further computation we have chosen the sim- 
plest form of the instability increment according to [l7j 
(weakly nonlinear Schrodinger equation with the value of 
the nonlinearity parameter e ~ 0.1 to 0.2): 



< / 



A0 

0Ak 



< 1. 



(2) 



Aiming to construct a stationary ("saturated") energy 
spectrum we regard cascading chain of the form 



200 = 01,1 + 02,1, 
01,1 = 02,1 + 02,2, 



0n,l 



3 n,2, 



Ei = piE , 
E 2 = p 2 Ei, 



(3) 



where < pj < 1, 0/ is the forcing frequency, Ej is the 
energy at the j-th step of cascade and pj is the part of 
the energy Ej-i transported from cascading mode Aj^i 
to cascading mode A,-. Each cascade step is regarded 
independently at discrete time moments. 

2.2.1. Assumptions. We assume further that 

(1) pj = p = const, i.e. cascade intensity p is constant 
for given excitation parameters; in this case p = (1 + 
VI — Eo)/2, where the energy of the initial excitation is 
normalized to 1 and Eq is the energy of its stationary 
amplitude Aq. 

(2) at each cascade step j, one cascading mode with 
frequency 0j,i is generated, with maximal increment of 
instability 

|(A0) m | 

imax,n . , ' KV 

(this corresponds to the idea of Phillips who suggested 
that in the saturated range the spectral density is sat- 
urated at a level determined by wave breaking). Modes 



0j 2 have smaller increments of instability and are fluxless 
("frozen", [loj]). 

2.2.2. The chain equation. Cascading chain consists of 
modes with frequencies 0i,i, 02, i, 03,1, it can be shown 
that the distance between two cascading modes |(A0)„| 
decreases slowly with growing n. It follows from (1) that 

E n =p n E A n+1 = y^A n , (5) 

The second assumption (2) can be rewritten as 

0n+l = 0n ± 0nA n k n , (6) 

and combination of A n+ i = ^/pA n and (O yields 

y/pA n = A(0 n ± n A n k n ) (7) 

which is further on called the chain equation ("+" and 
"-" stand for direct and inverse cascade, accordingly). 

2.2.3. Energy spectra. Taking Taylor expansion of 
RHS of the chain equation yields 

00 ^(«) 

ifpAn = A(0 n ± n A n k„) = ^ "HT (±0nA n k n ) S 

s=0 S ' 

= A n ± A n n A n k n + ^A^(±0 n A n k n ) 2 + ... (8) 
Taking two first RHS terms from ([5]) we get 

Vp- 1 



y/pA n — A n zb A n n A n k n A n — d=- 



0nkn 



(9) 



A{0) = ±yp -l)f^j- + const(0 o , A ) (10) 

Substitution of some specific dispersion relation into Q 
yields dependence A — A(0) and the form of energy spec- 
trum, E - A 2 . E.g. if = k a , a > 0, (HE) yields 



A{0) 



2 + a- 



?(0-( 2+ ^ 1 »-0o (2+Q " ) ) + A o , 



E(0) - 0~ v with 2 + a^ 1 < v < 2(2 + a" 1 ). 

This is in accordance with the laboratory results reported 
e.g. in 0,0], for ~ k 1/2 . 

2.2.4- Direction and termination of cascades. Condi- 
tions for cascade's termination and for formation of the 
direct and inverse cascades can be easily obtained from 
the form of the increment Namely, cascade termi- 
nates at some step N if is violated: 

0tv+i-0at = O or \0 N+ i - N \ > N A N k N . (11) 

Direct and inverse cascades occur if 

n+ i - 0„ > or 0„ - 0„ + i > (12) 

accordingly. Depending on the value of excitation pa- 
rameters, N < oo or N = oo, [19j j . 
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FIG. 2: Color online. Classical model of WT (on the left) and D-model of WT (on the right). 



A simplified form of cascade with |(A0)„| = const, 
20], though useful for interpreting some experimental 
results, yields cascade's termination at the final step and 
does not give a lead for possible transition to continuous 
spectrum. 

3. Classical model versus D-model. 

Graphical presentation of the classical model and D- 
model of WT is shown in Figj2] (in the left and right 
panels correspondingly). One of crucial assumptions of 
the classical model is, that forcing and dissipation are 
wide apart in the k-space, while the energy transport de- 
scribed by stationary spectra takes place within inertial 
interval (fcijfcg). This assumption is very difficult to re- 
produce either in laboratory experiments or numerical 
simulations. "Finite-size" effects are supposed to appear 
in the range of wavelengths k < k\ and are not described 
by the classical model, while in the range k > ki energy 
transport is supposed to be suppressed by dissipation. 
This phenomenology does not depend on the details of 
the initial forcing. 

On the other hand, in the D-model dissipation might 
be switched on or off at any cascade step by changing the 
magnitude of cascade intensity p, i.e. no additional as- 
sumptions are needed about relative positions of forcing 
and dissipation in the k-space; the existence of inertial 
interval is not important either and initial excitation 
may take place all over the k-space. In this model, three 
types of modes are singled out; energy transport in wave 
system depends on the type of the initially excited mode. 

3.1. Resonant and quasi-resonant modes (shown as red 
vertical solid T-shaped lines having a spring-like part); 
being excited, they may yield chaotic or recurrent pat- 
terns; no cascade shall be generated (e.g. In this 
case time synchronization condition (for frequencies) is 
satisfied exactly or with a small frequency mismatch 
while space synchronization condition (for wave vectors) 
is satisfied exactly: 

01+02 = 03 + 0, h + k 2 = k 3 , < fi < 1 or 

01 + 02 = 03 + 04 + h + k 2 = k 3 + fc 4 , 

and phases are coherent. 



3.2. Cascading modes (shown as small dashes compos- 
ing black dashed curves) ; being excited, they form a cas- 
cading chain with energy discrete spectrum ~ A 2 with A 
computed from (|10[) : for some specific initial magnitudes 
Aq and 0o ) energy spectrum might turn into continuous 
(shown as black solid "tail" of the curve) though not nec- 
essarily a classical spectrum. Choice of the parameters of 
initial excitation also defines whether excitation of a cas- 
cading mode yields direct or inverse cascade (e.g. [22J). 
In this case the time synchronization condition ([T]) is sat- 
isfied exactly, while the space synchronization condition 
is violated: 

0i + 02 = 203, ki + k 2 = 2fc 3 + O, < 6 < 1, 
and phases are are not coherent. 

3.3. Frozen modes (shown as blue vertical solid 
T-shaped lines): these are fluxless modes which do not 
take part in energy transport over the k-space, they 
may occur for some choice of initial state (e.g. [ToL l2lj|). 
For frozen modes both time and space synchronization 
conditions are violated. This means in particular that 
fluxless modes 01,2, 02,2, 03,2, ■•• have non-coherent phases 
and provide a necessary prerequisite in the statistical 
theory for a possible further development of a classical 
spectra. If at some cascade step a (quasi-) resonant 
mode is generated (shown as a black circle, FigHJ on the 
left) a cascading chain may be broken by appearance of 
an intermittency. 

It is important to realize that increment of instabil- 
ity is a parameter characterizing interplay of dispersion 
and nonlinearity in a wave turbulent system. Accord- 
ingly, two conditions of cascade termination (ITT|) have 
different physical consequences: 0at+i — 0n = yields 
1 = and manifests transition to a linear regime, while 
|0jv+i — 0n\ > 0NAjykjy corresponds to the case I > 1 
and growing nonlinearity. For instance, it has been shown 
in [19( for surface water waves that starting with a small 
nonlinearity Eq = Aq0q = 0.1, initial amplitude looses 
an order in the magnitude after 31 steps of the cascad- 
ing chains while the steepness of a wave packet becomes 
0.288. This means that another form of the increment 
should be used further on, e.g. 

I = \A \/(0Ak-^0 2 A 2 ky (13) 
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obtained in [18| for weakly nonlinear Schrodinger equa- 
tion with parameter of nonlinearity e ~ 0.25 to 0.4. This 
means also that the condition |0jv+i — 0n\ > 0n^-n^n 
allows to compute a critical frequency O rit (shown as a 
(black) bold vertical line at Fig., on the right) at which 
nonlinearity becomes mode important that dispersion. 
The energy spectra E{0), > cr u would differ substan- 
tially from the classical spectra. 

Another interesting phenomenon observed in numer- 
ous laboratory experiments is formation of a narrow 
zero-frequency sideband with non-zero energy (shown 
as vertical yellow rectangular). It should be checked 
whether this phenomenon might be attributed to termi- 
nation of an inverse cascade or originates in the (quasi) 
resonant interactions of short eigen-modes (appearance 
of an intermittency) . 

4. Conclusions. 

In this Letter our main goal was to understand natural 
phenomena revealed by laboratory observations of WT 
and to find the physical mechanism that underlies it. We 
have developed a simple model and analyzed it to see if 
it captures critical aspects of the phenomenon. 

We have demonstrated that D-model of WT allows to 
reproduce a number of discrete effects including gener- 
ation of a set of cascading modes in a wave turbulent 
system. We have also deduced the corresponding energy 
spectrum depending on the form of dispersion relation 
= 0(k), the form of the instability increment and pa- 
rameters of the initial excitation. Conditions for forma- 
tion of a direct or an inverse cascade have been written 
out explicitly. It has been shown that, depending on ini- 



tial state, a cascade might yield decrease or growth of the 
nonlinearity. 

A few possible ramifications of the model are accounted 
for below (the list is not exhaustive): 

• to change the form of the instability increment ([2]), e.g. 
to (fTUl) ; a number of examples is given in [l6[ ; 

• to regard cascade intensity as a function of wavelengths, 
p = p(k) 7^ const, then ODE similar to © is also solv- 
able in quadratures: 

A( ) = ± [ ( ^- 1)rf0 + const( 0o , Ac). (14) 
J 0k 

For instance, power function p(kj) = fcj may model 
dissipation; piecewise constant function p — pj , Vfe G 
(fej— l , kj) may model wind-generated surface water waves 
at different fetches (Jtj-i,kj), etc. 

• to regard the hierarchy of D-models with refined form of 
turbulent energy spectra obtained by cutting the Taylor 
expansion (J5J> at 3, 4 and so on terms instead of the 
simplest equation ©. It would be a challenge to find 
the general solution of the infinite order ODE given by 
([8]) and to investigate bounds of its physical applicability 
(it should be verified that energy flux corresponding to 
chosen ODE is not linear). 
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